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, Abstract. A review is given on the theoretical studies of charge correlations in 0- 

I— 1| (BEDT-TTF)2A'. Various studies show that within a purely electronic model on the 

^ ' 0-type lattice with the on-site U and the nearest neighbor Vp and Vc interactions, the 

, diagonal stripe, c-axis three-fold, and the vertical stripe charge correlations are favored 

' in the regime Vp < Vc, Vp ^ Vc^ and Vp > Vc, respectively. In the realistic parameter 

regime of Vc, there is a competition between c-axis three fold state and the 

diagonal stripe state. Since these are different from the experimentally observed a-axis 
three fold and the horizontal stripe charge correlations, additional effects have to be 
. included in order to understand the experiments. The electron-lattice coupling, which 

{Zj [ tends to distort the lattice into the fl^-type, is found to favor the horizontal stripe state, 

suggesting that the occurrence of this stripe ordering in the actual materials may not 
be of purely electronic origin. On the other hand, distant electron-electron interactions 
^ ' have to be considered in order to understand the a-axis three fold correlation, whose 

^ . wave vector is close to the nesting vector of the Fermi surface. These studies seem 

^sO ! to suggest that the minimal model to understand the charge correlation in 6'-(BEDT- 

' TTF)2Ar may be more complicated than expected. Future problems regarding the 

, competition between different types of charge correlations are discussed. 

(N ■ 

o ; 
o\ . 

O ■ random phase approximation, Fermi surface 

> : 

^ ■ 1. Introduction 

In recent years, charge correlations, i.e., orderings and fluctuations have been one of 
the central issues in the physics of organic materials [H El |3]. 6'-(BEDT-TTF)2X[4] is 
one of the most interesting families of organic compounds that exhibit short- or long- 
range charge correlations. The crystal structure of 0-(BEDT-TTF)2X, where BEDT- 
TTF = bis(ethylenedithio)tetrathiafulvalene, consists of a stack of alternating layers 
of BEDT-TTF and the anion X. In the BEDT-TTF layers, the molecules form an 
anisotropic triangular lattice, and the anion X controls the angle between the BEDT- 
TTF molecules, which in turn determines the band structure. Mori and coworkers [5l [6] 
obtained the electronic phase diagram in the angle-temperature space as shown in figure 
[H [7] which shows the wide variety of this series of materials, ranging from charge- 
ordered insulators such as X =RbZn(SCN)4 to the superconductor X =1^. Further 
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interest in this series of compounds has arisen by the observation of giant nonhnear 
transport for X =CsM'(SCN)4 (M' =Co, Zn), which causes the material to act as an 
organic thyristor [HI [9] . There, it has been pointed out that the coexistence of two types 
of short range charge correlations with modulation wave vectors {qa,Qb,(lc) = (f,^, |) 
and (0, k, |) (in units of the reciprocal lattice primitive vectors in the a, b, and c 
directions; the a and c directions are shown in figure [2]) plays an important role in 
this phenomenon. The (0, k, |) charge modulation is often referred to as horizontal 
stripe ordering since this wave vector corresponds to a stripe-type charge modulation, 
where the stripes run in the horizontal direction (figure [3]). The (|,/c, |) modulation 
corresponds to a 3 x 3 charge correlation with respect to the a-c unit cell, as shown in 
figurelH At high temperatures, a diffuse X-ray spot is observed at (|, k, |) [10], [11], while 
the (0, k, |) structure develops as the temperature is lowered, and the system becomes 
more insulating (figure [S]). When an electric field is applied, the (0, /c, |) horizontal 
stripe ordering is degraded, resulting in a recovery of the metallic behavior and thus the 
nonlinear transport. Related to this phenomenon is the observation of lattice modulation 
with wave vector (|, k, 0.29 -|) for X =CsCo(SCN)4 at a pressure of 10 kbar, which has 
been attributed to a pressure induced 2kF charge density wave because the modulation 
wave vector coincides with the nesting vector of the Fermi surface [TT] . 

For X =RbM'(SCN)4 (M' =Co,Zn) also, two types of charge correlation are 
involved. At high temperatures in the metallic phase, diffuse X-ray spots are observed at 
(|, k, i) and (|, k, |), which suggests the presence of short-range 3x4 charge correlation 
(figure Hj) [121 US]- Some anomalies have also been observed in the nuclear magnetic 
resonance experiments in a similar temperature range [HJ [1^. At around 200 K, 
the system undergoes a metal-insulator transition, accompanied by a structural phase 
transition to the so-called 6*^ phase [16], in which the unit cell is doubled in the c 
direction [6] [121 (ISl O (HI [13 [HI [19]. In the 9d phase. X-ray diffraction measurements 
have revealed that long-range horizontal stripe charge ordering with the modulation 
wave vector (0, 0, ^) takes place [5l [6l [121 [13] . 

Theoretically, the origin of these charge correlations has been an issue of great 
interest [2]. A number of analytical and numerical studies have been performed using 
the extended Hubbard model (coupled with the lattice in some cases) described in the 
following section. In the present paper, we review theoretical studies on the charge 
correlations in ^-(BEDT-TTF)2X. 

2. Theoretical model 

In this section, we describe an electronic model of ^-(BEDT-TTF)2X, where each 
molecule is considered as a "site" in the tight binding approximation. In the lattice 
structure of ^-(BEDT-TTF)2X shown in figure [21 the direction of the molecules 
alternates along the a axis resulting in the rectangular a-c unit cell, but this alternation 
is irrelevant as far as the hopping integrals in the tight binding model are concerned. 
Therefore, we can consider a unit cell (x-y) which is half the size of the usual unit cell. 
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Figure 1. Mori's phase diagram for 0-(BEDT-TTF)2X. Reproduced with permission 
from ref. [7| . Copyright 2006 by the Physical Society of Japan. 



thereby obtaining the effective lattice structure shown in figure El The corresponding 
Brillouin zone becomes unfolded, and the relation between the unfolded and folded 
Brillouin zone is shown in figure [71 Here, the wave vectors are given in units of 
the reciprocal lattice primitive vectors, and the subscripts uf and / indicate whether 
the reciprocal vectors are those of the unfolded or folded Brillouin zone, respectively. 
The wave vector of the experimentally observed horizontal stripe charge ordering is 
(1/4, l/4)uj in the unfolded Brillouin zone. The wave vector (|, fc, | ~ i)/ at which the 
diffuse X-ray spots are observed in the high temperature regime lie on the diagonal line 
labeled as ka = |, which satisfies kx — ky = | in the unfolded Brillouin zone. Since these 
wave vectors correspond to charge correlation with threefold periodicity in the a-axis 
direction (figure Hj), we will call these the "a-axis threefold" charge correlation. 

Important progress in understanding the origin of the charge correlations in 6- 



Theoretical Aspects of Charge Correlations in 9-(BEDT-TTF)2X 



4 




Figure 2. Lattice structure of the cation layer of 6'-(BEDT-TTF)2X. tp and tc are 
the nearest neighbor hopping integrals. {tp,tc) — (106, —5), (137, —49), and (99, —33) 
[meV] for CsCo, CsCo under pressure (10 kbar), and RbCo salts, respectively, according 
to the extended Hiickel estimation [TTj. Vp and Vc are the nearest-neighbor interactions, 
while Va, Vq, and V2c are the next-nearest- neighbor interactions. The a-c unit cell is 
the usual unit cell, while we can use the x-y unit cell to unfold the Brillouin zone. 



(BEDT-TTF)2X was made by Seo, who included the nearest neighbor off-site repulsions 
Vp and Vc in the model Hamiltonian (see figure [2])p6| in addition to the on-site 
interaction U, which had already been considered in the seminal studies of organic 
conductors by Kino and Fukuyama [20] and Seo and Fukuyama [2T]. Mori estimated 
the interaction parameters by directly calculating the Coulomb interaction using the 
highest occupied molecular orbitals of BEDT-TTF, neglecting the screening effect |22j . 
This estimation shows that Vp ^ Vc {Vp is slightly smaller than Vc), while interactions 
Va, Vp, and V2c are also of significant size compared with Vp and Vc- 

A number of theoretical studies on the charge correlations in ^- (BEDT-TTF) 2X 
have been performed using the extended Hubbard model on the lattice shown in figure 
|6] [23] , where the Hamiltonian is given in the form 

UjC^Cj^ + Uj2'^i]'f^il + J2 ^ij^i^j^ (1) 



H 



E 

<i,j>,a 



where c]^ creates an electron with spin a =|, J, at site i, and = cj^Cio-. The band 
filling, i.e., the average number of electrons per site (molecule), is fixed at n = 1.5 in 
accord with the actual material. In figure O the letters p, c, a, ■ ■ ■ denote the relative 
positions with respect to the red site at the center, and we use these letters to specify 
the range of the hopping integrals and the off-site repulsive interactions. For the kinetic 
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Figure 3. Horizontal stripe, vertical stripe, diagonal stripe, and c-axis threefold charge 
correlations. The green (or dark) molecules are the charge-rich ones. 
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Figure 4. a-axis threefold (3x3 and 3x4) charge correlations. The dark molecules 
are the charge-rich ones. 



energy terms, we consider hoppings t^ and t^-, where tp (~ 0.1 eV) is taken as the unit of 
energy. X =CsM'(SCN)4 (M' =Co, Zn) has small values of t^, while X =RbM'(SCN)4 
[M' =Co, Zn) and X =CsCo(SCN)4 under pressure have relatively large values of 
~ -0.3 to -0.4 [n]. 

3. Stripe-type charge correlations 

Three types of stripe-type charge correlations, vertical, horizontal, and diagonal (figure 
|3]), were proposed by Seo [16]. The wave vectors of the three patterns correspond 
to (l/2,l/2)„; = (0,1)/, (l/4,l/4)„; = (0,1/2);, and (l/2,0)„y = (1/2,1/2);, 
respectively. Intuitively, the horizontal and diagonal stripe states are favored when Vp 
is small, while the vertical stripe state is favored when Vc is small, as can be understood 
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Figure 5. Upper panel: X-ray scattering intensity scanned along a diffuse sheet, 
expressed in units of the reciprocal lattice primitive vectors as (3 1 1/2) + (2^ £,) 
with —0.09 < ^ < 0.48, where the horizontal axis in the figure is h = —3 + 2^. qi 
is the wave vector corresponding to the a-axis threefold (3 x 3) correlation, while (72 
corresponds to the horizontal stripe correlation. Reproduced with permission from ref. 
[TT] . Copyright 1999 by the Physical Society of Japan. Lower panel: line (the arrow) 
along which the X-ray scan is taken (i.e., the diffuse sheet) shown in the Brillouin zone 
adopted in theoretical studies (see figure [7|). 



from figure [3l In ref. [16] , the energies of various types of charge correlation patterns 
were calculated within the mean field approximation. The calculation shows that the 
diagonal stripe is more stable than the horizontal and vertical stripe states when Vp is 
small, while the vertical stripe state is stable for small Vc- Exact diagonalization studies 
on these stripe orders have also been performed [231 125] . 

Seo further showed that the horizontal stripe state is much stabler in the 9d phase 
of the crystal structure than in the 9 phase, where the hopping integrals are modulated 
as shown in figure M [26] In this 6d lattice structure, one-dimensional chains of molecules 
connected by the large hopping p4 are weakly coupled by other small hoppings. In 
Seo's calculation, the competition between the diagonal and horizontal stripe states was 
subtle, but the resulting effective spin Hamiltonian in the horizontal stripe phase is a 
quasi-one-dimensional Heisenberg chain along the pA hopping integral, which explains 
the Bonner-Fisher type temperature dependence of the magnetic susceptibility [5l E]. 
From this viewpoint, Seo concluded that the horizontal stripe ordering is the most 
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Figure 6. Effective iattice structure using the x-y unit cell. We specify the hoppings 
and the interactions by c, p, a, q, and 2c, which specify the positions relative to the 
center (red) site. 



probable state in the insulating phase of X =RbM'(SCN)4. This was later confirmed 
by X-ray experiments [I2l |T3] as mentioned in the introduction. 



4. Competition between stripe and threefold charge correlations 
4.1. c-axis threefold state 

After the seminal study by Seo Mori performed an analysis in the static limit (where 
the band structure is neglected), and pointed out the possibility of the nonstripe-type 
charge correlation shown in figure [3] [22] , which has threefold periodicity in the c-axis 
direction. We will call this the "c-axis threefold" charge correlation. The corresponding 
wave vector is (1/3, l/3)u/ = (0, 2/3) /. In real space, the c-axis threefold state is favored 
when ~ and Va = Vg = V2c = because this configuration minimizes the energy 
cost of Vp and as can be seen from figures [2] and [31 

Kaneko and Ogata extended Seo's mean field study to take into account the 
possibility of this c-axis threefold state [27]. They found that this state has a lower 



energy than the vertical, diagonal, and horizontal stripe states when Vp ^ Vc and 
Va = Vq = V2c = 0. Hotta and coworkers used exact diagonalization for a spinless 
half-filled model, which can be considered as an effective model of ^-type compounds 
in the large on-site U limit, and they showed the presence of c-axis threefold charge 
correlation for K ^ K [28[l29]. 

A variational Monte Carlo study by Watanabe and Ogata [30] , and a later density 
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Figure 7. Upper panel: relation between the folded and unfolded Brillouin zone. 
Modulation wave vectors corresponding to various correlations are shown. The wave 
vectors are presented in units of the reciprocal lattice primitive vectors. The red line 
shows the Fermi surface for = —0.3. Lower panel: Brillouin zone in the upper panel 
squeezed in the = —ky direction, and then placed repeatedly. The dashed lines 
represent the diffuse sheets observed at room temperature and below. 



matrix renormalization group study by Nishimoto et al. [3T] revealed the phase diagram 
of the model in Vp-Vc space with = and U/tp = 10. In the variational Monte 
Carlo study, the energy of each ordered state was calculated and compared, while 
in the density matrix renormalization group study, the real-space charge modulation 
was directly observed. As shown in figure [9l the two approaches give similar results, 
where the c-axis threefold state appears in a regime with Vp ^ Vc between the diagonal 
(Vp < Vc) and vertical stripe {Vp > Vc) regimes. 

A random phase approximation (RPA) analysis gives a clear view, in the reciprocal 
space, as to why the c-axis three fold state is favored in the Vp Vc regime, and also 
why the diagonal stripe is more favored than the horizontal stripe [32]. Within the RPA 
[33t [Ml [35] , the charge susceptibility Xc is given as 

r . _ X0(g) (r)^ 

l + {U + 2V{q))xo{qy ^ ^ 
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Figure 8. 9d lattice. The hopping integrals are cl — 1.5, c2 = 5.2, pi = 16.9, 
p2 = -6.5, p3 = 2.2, and p4 = -12.3 in units of 10"^ gy [12]. The green molecules 
are the charge rich ones in the horizontal stripe state. 



Here Xo is the bare susceptibility given by 

xo{q 



(3) 



p ^P~^p + q 
with ej^ being the band dispersion given as 

ej^ = 2tp[cos(A;^) + cos{ky)] + 2tcCos{kx + ky) (4) 

and /(ep) = l/(exp(ep — ij)/T) + 1) is the Fermi distribution function. When the 
nesting of the Fermi surface is good, Xoio) is maximized at the nesting vector. V{q) is 
the Fourier transform of the off-site repulsions, given as 

V{q) = 2Vp[cos{q^) + cos(gy)] 

+ 2Vc cos(ga; + Qy) + 2K cos(g^ - qy) 
+ 2Vq[cos{2qr, + Qy) + cos(g^ + 2qy)] 

+ 2V2cCOs{2q^ + 2qy) (5) 
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Figure 9. Ground-state phase diagram obtained by (a) variational Monte Carlo and 
(b) density matrix renormalization group for ic = 0, [/ = 10, and Va = Vq = V2c = 0. 
In (b) , the diameter of the circles is proportional to the hole density at the hole- rich 
sites, (a) Reproduced with permission from ref. [30j. Copyright 2006 by the Physical 
Society of Japan, (b) Reproduced with permission from ref. [31] . Copyright 2008 by 
the American Physical Society. 



In figure [ini the RPA charge susceptibihty for tc = —0.3, U = 3, Vc = 1.5, 
Vp = 1.3, and Va = Vg = = is shown. This choice of tc corresponds to the case of 
X =RbM'(SCN)4 (M'=Zn, Co) or X =CsCo(SCN)4 under pressure. One notices that 
there is a peak at (1/3, l/3)u/, which corresponds to (0,2/3)/ in the original Brillouin 
zone, namely, the modulation wave vector of the c-axis threefold correlation. This is 
due to a peak in —V{q) near q = (1/3, 1/3)^/ as shown in figure fTOl Namely, although 
Xo{q) is maximized around (2/3, l/4)j = (11/24, — 5/24)„j, the effect of —V{q) is strong 
in the denominator of the RPA charge susceptibility eq.(l2]), thereby minimizing the 
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denominator and maximizing Xc{q) (1/3; In addition to the c-axis threefold 

peak, a subdominant peak also appears around the wave vector close to (1/2, 0)„/, i.e., 
the wave vector of the diagonal stripe. This is due to a combination of two effects: (i) 
Xo{q) takes large values in the region extending from the nesting vector position toward 
(1/2, 0)u/, and (ii) the region where —V{q) takes large values tends to extend toward 
(1/2, 0)u/, particularly when Vp < Vc- A similar result is also obtained for smaller 
values of tc, although the subdominant peak near (1/2, 0)^/ becomes relatively weaker. 
Such a two-peak structure has also been found in the dynamical charge correlation 
function calculated using exact diagonalization, where peaks at low energies are found 
at (1/2,0),/ and (1/3, l/3)„/ [3l]. 

Udagawa and Motome [36] further investigated the temperature dependence of these 
peak structures in the charge susceptibility, and showed that the (1/3, 1/3)^/ peak is 
nearly temperature- independent, while the peak near (1/2, 0),/ increases in magnitude 
as the temperature is lowered. This is because the (1/3, 1/3),/ peak is mainly due to 
the nearest-neighbor electron-electron interaction, while the (1/2, 0)„/ peak is related to 
the Fermi surface nesting, which becomes more effective at low temperatures. From this 
temperature dependence, they claimed that these peaks correspond to the two diffuse 
X-ray spots observed in X =CsM'(SCN)4 (M' =Co, Zn) [TOl [11] (figure ED, although 
the wave vectors are different. 

4-2. Effect of electron- lattice couplings 

As seen in the previous section, a purely electronic model on the ^-type lattice fails to 
reproduce the horizontal stripe charge correlation observed experimentally. Following 
the line of Seo's study fT6], Kaneko and Ogata [27] performed a mean-field analysis on 
the 6d lattice structure, where competition between the c-axis threefold state and the 
horizontal stripe was found to take place. 

Tanaka and Yonemitsu took into account Peierls-type electron-lattice couplings, 
and performed a Hartree-Fock calculation [37|. They showed that the horizontal stripe 
state is stabilized by a self-consistently determined lattice distortion into the 6'^^ phase 
when Vp/Vc < 1. The stabilization of the horizontal stripe state by electron-lattice 
coupling has been further supported by an exact diagonalization study by Miyashita 
and Yonemitsu [38] . 

Udagawa and Motome included the Su-Schrieffer-Heeger type electron-phonon 
interaction, which also takes into account the effect of the lattice distortion, and 
performed a Hartree-Fock calculation [36]. They found that the horizontal stripe state 
is stable for sufficiently large electron-phonon couplings. 

Note that although these calculations succeed in reproducing the horizontal stripe 
charge correlation, the competing candidates considered therein do not include the a- 
axis threefold correlation observed experimentally. 
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Figure 10. Contour plots of the RPA result of the charge susceptibility Xc, the 
bare susceptibility xo: and the Fourier transform of the off-site interactions —V{q). 
tc = -0.3, [/ = 3, Vp = 1.3,14 = 1-5, Va = Vg = V2c = 0, and T = 0.05. The 
solid (dashed) arrow in Xc{<l) indicates the dominant (subdominant) peak at wave 
vectors close to (1/3, l/3)„/ ((1/2, Q)uf)- The solid arrow in —V{q) is the wave vector 

(i/3,i/3)„/ m- 



4-3. a-axis threefold charge correlation 

As seen in the previous subsections, models that consider only U, Vp, and Vc as electron- 
electron repulsions have a tendency toward the c-axis threefold state formation in the 
realistic parameter regime of ~ Vc- This state competes with the diagonal (horizontal) 
correlation in the absence (presence) of electron-lattice coupling [39]. Correspondence 
between the c-axis threefold charge correlation and the experimentally observed short 
range a-axis threefold charge correlation has often been discussed, but the wave vectors 
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Figure 11. Contour plots of the RPA result for the charge susceptibility Xc^ the 
bare susceptibility xo, and the Fourier transform of the off-site interactions —V{q). 

= -0.3, U = A,Vp = 1.8, 2.0, K = 0.4,^, = 0.7, V2c = 1.1, and T = 0.25. 
The arrow in xo(q) indicates the peak position, i.e., the nesting vector of the Fermi 
surface. The arrows in —V{q) indicate the area at which ~V{q) is broadly maximized 
|32|. 



of the two threefold correlations are different. 

Kuroki performed an RPA analysis on a model that considers the distant 
interactions V^, Vg, and in order to investigate the origin of the a-axis threefold 
correlation |32]. In figure [TT], the charge susceptibility for t^ = —0.3 is shown, where a 
peak is found near the a-axis threefold correlation wave vectors. Here only the result for 
certain values of Va, Vq, and V2c is shown, but in ref. [32] it was found that similar results 
are obtained within a certain range of Va, Vq, and V2c values that are smaller than ~ 
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and ~ Vc/2. This result can be understood as follows. In the presence of the distant 
interactions, —V{q) is now broadly maximized in a region that includes the nesting 
vector position, as shown in figure [TTl Therefore, the denominator in eq.([2]) is now 
minimized at a position close to the nesting vector, resulting in a charge susceptibility 
peak position that nearly coincides with the diffuse X-ray spot positions. As mentioned 
in the introduction, the relevance of the Fermi surface nesting to the a-axis threefold 
correlation had already been pointed out in ref. [H]. Thus, the a-axis threefold charge 
correlation, within this scenario, should be considered as being conceptually different 
from the c-axis threefold charge ordering in that it is not related to the Fermi surface. 

In figure [121 the result for t^ = —0.05 with the other parameters having the same 
values as those in figure [TT] is shown. This choice of tc corresponds to X =CsM'(SCN)4 at 
ambient pressure [11]. In this case, the peak position of the charge susceptibility moves 



Theoretical Aspects of Charge Correlations in 6-(BEDT-TTF)2X 16 

4 = -03t, (Rb) 

71 



~-n 271 371 



= -O.OStp (Cs) 

7t 



-71 JC 271 37C 

Figure 13. Comparison of tlie Fermi surface nesting between the cases of tc = —0.3 
and -0.05. 



0.1 



f4 0.01 



0.001 

-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 

^ tA 

smaller deviation between -V(q) and max. positions 
better Fermi surface nesting 
larger density of states at £f 
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X =13, the value of tp is taken from ref. [B], while that of tc is from ref. [1T| . 
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toward the unfolded Brillouin zone edge, i.e., closer to (|, — = (|)|)/; compared 
with the case when t^. = —0.3. This is because the nesting vector is closer to (|, |)/, as 
can be seen in the figure of xo, which in turn can be understood from the Fermi surface 
shown in figure [T3l This tendency appears to be consistent with the experimental finding 
that diffuse X-ray spots at high temperatures are observed at Qc = ^ iia. X =CsM'(SCN)4 
at ambient pressure, while they are at | < gc < | in X =CsM'(SCN)4 under pressure 
and in X =RbM'(SCN)4. 

In figure [T^ the temperature T^-axis at which the RPA charge susceptibility diverges 
at a certain q is shown as a function of tc [32j. T^-axis decreases with increasing tc 
(from tc < to tc > 0), and sharply drops when tc > 0. In fact, the increase in tc 
induces several effects. First, the nesting of the Fermi surface is degraded (see figure 
[T3l) . second, the nesting vector deviates from the region where —V{q) is maximized, 
and third, the density of states near the Fermi level becomes smaller. These factors 
together deteriorate the charge ordering. The decrease in the density of states with 
increasing tc has also been pointed out in ref. [ID]. Figure [TH is reminiscent of the phase 
diagram of the ^-(BEDT-TTF) family obtained by Mori et al. [6], where charge ordering 
takes place for compounds having large negative tc such as X =RbM'(SCN)4, while it 
does not occur for compounds having small or positive tc [HI SSI US] such as X =1^. 
However, note that in figure [TH the charge susceptibility diverges at a modulation wave 
vector close to the nesting vector of the Fermi surface (around (|, | — |)/ for tc < 0), 
while the charge ordering temperature in the experimental phase diagram is that for 
the (0, fc, |)/ ordering, i.e., the horizontal stripe state. The resemblance between the 
two figures suggests a relationship between the horizontal stripe and the a-axis threefold 
correlations. This may also be consistent with an experimental analysis under uniaxial 
pressure, which concluded that the value of tc is the key factor that dominates the 
occurrence of horizontal charge ordering [Uj. A possible relation between the a-axis 
threefold and the horizontal stripe orderings will be mentioned in the next subsection. 

4.4- Future problems regarding the competition 

Since the distant electron-electron interactions Va, Vg, and V2c are found to be necessary 
to reproduce the experimentally observed a-axis threefold charge correlations, it is 
reasonable to assume that the related experimental observations have to be understood 
within such a model. 

One problem to be resolved is the temperature dependences of the two diffuse X-ray 
peaks observed in X =CsCo(SCN)4 [TT]. Since the peak in the charge susceptibility 
partially originates from the Fermi surface nesting in the scenario described in 
the previous subsection, the maximum value of the peak increases with decreasing 
temperature, which is inconsistent with the experimental observation that the X-ray 
peak corresponding to the a-axis threefold charge fluctuation is nearly temperature- 
independent below 90K [llj. Regarding this problem, it is possible that the anions, 
which are not taken into account in previously studied models, may play some role. 
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Figure 15. xo{q) for tc = -0.05 at T = 0.15 (- 150 K) (left) and at T = 0.02 (- 20 
K) (right). 

Namely, in the X-ray scattering experiment, diffuse sheets (or planes) along the 2a + c 
and the 2a — c directions are observed at room temperature below, which is considered 
to be due to the effect of anions (the diffuse sheets are shown in the lower panels of figure 
[5] and figure [7]) plJ]. Actually, the wave vectors (|, k, |) (a-axis three fold) and (0, k, i) 
(horizontal stripe) are on these diffuse sheets, and they "grow" from these sheets as 
the temperature is lowered, suggesting the relevance of the anions. Kuroki considered 
the effect of these diffuse sheets phenomenologically by adding a term that enhances the 
charge susceptibility along the hne 2a+c to the denominator of eq.([2]), and he found that 
the maximum value of the charge susceptibility saturates at low temperatures [15]. This 
is because the wave vector at which Xo{q) takes its maximum gradually moves away from 
(2/3, 1/3)/ toward (3/4, l/4)j as the temperature is lowered as shown in figure [15] (note 
that the nesting is less good for the case of X =CsCo(SCN)4), so that the cooperation 
between the diffuse sheets and the Fermi surface nesting becomes degraded. 

Another problem is the origin of the horizontal stripe charge correlation. As 
mentioned in the previous sections, in the absence of the distant interactions Va, Vq, 
and V2c, the horizontal stripe state is shown to be favored by electron-lattice coupling 
that distorts the lattice structure into the 6'd-type. Nevertheless, it remains a future 
problem as to whether this conclusion holds in the presence of the distant interactions. 
Moreover, even if the origin of the horizontal stripe in the Rb salt is the lattice distortion, 
the development of the horizontal stripe correlation at low temperatures in the Cs 
salt, where such a distortion has not been observed experimentally, remains puzzling. 
Related to this point, Kuroki pointed out the possibility that the horizontal stripe 
charge correlation arises from the a-axis threefold correlation [32]. For the case of the 
Cs salt, although the nesting vector is close to (2/3,1/3)/ at high temperatures, it 
moves toward (3/4, 1/4)/ (note that we are adopting the folded Brillouin zone scheme 
here) upon lowering the temperature as mentioned above. Thus, at low temperatures, 
we have two different nesting vectors close to (3/4,1/4)/ and (—3/4,1/4)/. Their 
coexistence in the momentum space might result in a modulation having the wave 
vector (0, 1/2)/ = (3/4, 1/4) + (—3/4, 1/4), which corresponds to that of the horizontal 
stripe. This may be considered as a kind of "4^^?" correlation in that the wave vectors 
(3/4, 1/4)/ and (—3/4, 1/4)/ are close to the nesting vectors of the Fermi surface, namely, 
the "2/cp'" wave vectors at low temperatures. Investigation of the occurrence of the 
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horizontal charge correlation due to the interplay between such a correlation and 

the electron-lattice couplings (including the interaction with anions) will serve as an 
interesting future area of research. 

There are several studies in which the pairing symmetry of superconductivity in 
X =13 has been discussed [301 [321 [351 US] . Since this superconductivity lies in the vicinity 
of the charge-ordering phase in Mori's electronic phase diagram (figure [T|), one may 
expect the occurrence of unconventional superconductivity due to charge fluctuations. In 
such a case, the form of the gap function strongly depends on the wave vector dependence 
of the charge fluctuation. For instance, in a model that considers electron-electron 
interactions up to nearest neighbors, the coexistence of spin and charge fluctuations 
results in competition between spin-singlet cij^y-wave-like pairing and triplet /-wave-like 
pairing [30l [35]. On the other hand, in the model that takes into account the next 
nearest neighbor interaction, the gap function can be more complicated [32]. Since 
it can be expected that the many body interactions do not depend strongly on the 
anions (although the band structure does), one must discuss the pairing symmetry for 
X =13 using a Hamiltonian with appropriate interactions that correctly reproduces the 
competition between the a-axis threefold and the horizontal stripe charge correlations in 
the case of Rb and Cs salts. In this sense, the superconducting pairing symmetry should 
be theoretically determined after the origin of the charge correlations is completely 
understood. 

5. Conclusion 

In the present paper, we have reviewed theoretical studies on the charge correlations 
in ^^-(BEDT-TTF)2MM'(SCN)4. It now seems clear that within the purely electronic 
model on the 6'-type lattice with the on-site U and the nearest neighbor Vp and Vc 
interactions, the diagonal stripe, c-axis threefold, and vertical stripe charge correlations 
are favored in the regimes Vp < Vc, Vp Vc, and Vp > Vc, respectively. In the realistic 
parameter regime of ~ Vc, there is competition between the c-axis threefold state and 
the diagonal stripe state. Since these are different from the experimentally observed a- 
axis threefold and the horizontal stripe charge correlations, additional effects have to be 
taken into account to understand the experimental results. The electron-lattice coupling, 
which tends to distort the lattice into the ^^-type, is found to favor the horizontal stripe 
state, suggesting that the occurrence of this stripe state in actual materials is not of 
purely electronic origin. 

On the other hand, the inclusion of the electron-lattice coupling does not explain 
the presence of the short-range a-axis threefold charge correlation, and we have to 
add distant (next-nearest-neighbor) electron-electron repulsions to understand this 
phenomenon. Related to this point, a transport experiment has indicated that the 
repulsive interaction in ^-(BEDT-TTF)2MZn(SCN)4 (M =Cs,Rb) is long-range [37]. In 
fact, the presence of such a distant interaction has also been pointed out theoretically 
[m [321 [51] for quasi-one-dimensional organic compounds (TMTSF)2X in the context 
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of coexisting 2kp spin and 2/cf charge density waves [521 E3] and spin triplet pairing 
superconductivity [501 [511 [SH |55l [56], |57l [58] . Quite recently, Tahara et al. estimated 
the electron-electron interaction in fi;-(BEDT-TTF)2X from a first principles calculation, 
and also obtained a long ranged interaction [59] . Thus, it is possible that the long-range 
nature of the electron-electron interaction may be common in these organic materials. 

So far there have been few theoretical studies that take into account the distant 
electron-electron interactions; thus, it is not clear whether competition between the a- 
axis threefold and horizontal stripe states takes place in the presence of electron-lattice 
couplings in such models. At the present stage, one cannot completely rule out the 
possibility that strong correlation effects that are not taken into account in RPA or 
mean-field analysis may result in the horizontal stripe state in the presence of distant 
interactions. Also, as mentioned in section l474l it is possible that the anions, which have 
not been taken into account in theoretical models so far, play some role in the occurrence 
of these charge correlations. The minimal theoretical model may be more complicated 
than expected, and more theoretical study must be carried out to fully understand the 
physics of charge correlations in 6'-(BEDT-TTF)2X. 
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